Design of a technique to measure the density of ultracold atoms in a short-period optical lattice in 

three dimensions with single atom sensitivity 
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A measurement technique is described which has the potential to map the atomic site occupancies of ultracold 
y—i • atoms in a short-period three-dimensional optical lattice. The method uses accordion and pinning lattices, 

t-H ' together with polarization gradient cooling and fluorescence detection, to measure the positions of individual 

, atoms within the sample in three dimensions at a resolution of around half the atomic resonant wavelength. The 

s ! ■ method measures the site occupancy, rather than the parity of the site occupancy, of atoms in the lattice. It is 

5_l ' expected that such measurements hold significant potential for the study of ultracold quantum dynamics. 

In recent years there has been a surge of interest in the exploration of strongly interacting quantum systems of ultracold atoms 
, in optical lattices. In addition to long coherence times, a high degree of control over the quantum Hamiltonian is possible, 
' making such systems ideal for studies of many -body entanglement (2]. The direct measurement of the number of atoms at 
, each site of the optical lattice provides a great deal of information about the system, particularly in terms of measurements of 
& ' local correlation and excitation properties. The measurement of the positions of individual atoms within a strongly interacting 
quantum system in an optical lattice has been demonstrated very recently in two experiments JHI^t]; it may be anticipated that 
I ■ the excellent results from these experiments will stimulate further research in this direction. 

The present article outlines and models a measurement technique which potentially offers several advantages over the mea- 
' surement techniques employed in the recent experiments JUS]. The aim of the technique is to measure the number of atoms at 
each lattice site, for a dense three-dimensional sample of ultracold atoms held in a three-dimensional optical lattice formed by 
counter-propagating laser beams, for all lattice sites within a certain small volume (of order 10 x 10 x 10 lattice sites to a side). 

The technique which will be outlined in this article, in common with the recent experiments JUIS], uses fluorescence detection 
to measure the atomic density distribution. There are four major issues which must be overcome when using fluorescence 
I . detection to probe dense neutral atomic samples in an optical lattice on an atom-by-atom level. Firstly, each atom must be 
' O ' kept (to a high probability) within a small volume during the photon collection time in order to extract a strong and consistent 
spatially-localized signal. A deep optical lattice can be used to confine the atoms during the measurement, with the strong photon 
recoil heating compensated with polarization gradient cooling. This technique has been used in the recent experiments Qlit], 
and is explored in detail in a separate article by the author J4|]. 

Secondly, the spatial resolution of the image should be high enough to resolve individual lattice sites. Strongly interacting 

■ dynamics occur preferentially in small-period lattices and at high atomic density, leading to substantial difficulties in resolving 
0\ | the sites; nevertheless this capability has been successfully demonstrated in the two recent experiments J2, 0] using very high 

■ numerical aperture lens systems. 
Thirdly, exposure to intense resonant light during the measurement causes strong light-assisted collisional losses. The energy 

released in such collisions is such that the atoms are immediately ejected from the confining potential. These losses occur 
pairwise; isolated atoms are immune. The atoms are lost in a very small fraction of the measurement time, leaving no identifiable 
fluorescence signal. The current generation of resolved-atom experiments |0 0] measure the parity of the lattice site 
occupation numbers rather than the true number distribution. 

A fourth issue occurs when making measurements of a three-dimensional sample; such a measurement typically involves the 
tomographic reconstruction of the spatial distribution from a sequence of images at different depths. The measurement technique 
should be capable of distinguishing the signal from atoms at different depths within the sample; scattering from out-of-focus 
atoms should not overwhelm the signal from in-focus atoms. 

This article outlines a technique which has the potential to resolve the latter three issues. 

The accordion lattice is an optical lattice with a dynamically and smoothly variable spatial period J3-01- In Ref. JUl it is 
outlined how the lattice period can be controlled electro-optically using an acousto-optical deflector. The final optical element 
in the accordion lattice apparatus takes two beams, both parallel to the optical axis, directing them onto the atoms; the accordion 
lattice periodicity depends on the off-axis displacement s of the beams. When the final optical element is a lens with numerical 
aperture NA (Fig. QI a )X me minimum lattice period is limited to d m i n = A/(2 NA). Better performance is achieved using a 
parabolic mirror in place of the lens (Fig. Q2b)); this reduces the minimum lattice period to the optimal d^ n = A/2. A three- 
dimensional 'parabolic' accordion lattice can be formed using three independently controlled frequency-offset one-dimensional 
accordion lattices together with an appropriately constructed mirror element. One conceptual three-dimensional parabolic lattice 
mirror design is given in Fig. [2] 

The maximum spatial period of the accordion lattice formed using a parabolic mirror of focal length / is d mix = fA/{2s m i n ). 



(a) 



(b) 



FIG. 1: The basic geometries for the accordion lattice; only the last optical element is shown. The laser beams enter the last 
element parallel to the optical axis with an equal and opposite off-axis displacement which can be electro-optically controlled 

(see Jit] for details). The lattice spatial period in either configuration is d — A/ (2 sin 0). Here 6 is the angle between each laser 
beam and the optical axis at the position of the atoms; the maximum angle 6 max possible for the lens of Fig. (a) is given by the 

numerical aperture of the lens as sin ma x= NA. 




FIG. 2: (Color online) A mirror design for a three-dimensional parabolic counter-propagating accordion lattice. The underlying 
concept is to use a cylindrically-symmetric parabolic reflector (which has a central slice removed) to form a two-dimensional 
accordion lattice; a separate slim parabolic reflector is used to form the accordion lattice along the remaining dimension. Each 
pair of beams is generated separately. Unused parts of the parabolic reflector can be removed to improve the optical access. It is 
anticipated that the atoms would be observed using an objective lens (not shown) located in a diagonal position. The displayed 

design is not unique; variants are possible. 

This is determined by the closest distance s m [ n the incident beams can come to the optical axis (Fig. 02b)) without encountering 
a mechanical obstruction (e.g. see Fig. |2j or overlapping with the atoms prior to reflection from the mirror. 

In an experiment, the coherent quantum dynamics would occur while the lattice is in the minimum spatial period configuration, 
as both tunneling and interaction dynamics are strongest at small lattice periods. After the coherent dynamics have occurred, the 
lattice depth is increased to shut off interwell dynamics, pinning the atoms in place. The spatial period of the accordion lattice 
would then be smoothly increased for resolved-site measurements of the atomic occupation numbers. 

The parabolic mirror (Fig. |TJb)) both deflects and focuses the incident laser beams; collimated incident beams will be brought 
to a focus at the atoms. To compensate for this, the incident beams should have a focus before the mirror (on the front focal sur- 
face), so they are collimated at the atoms. However, the front focal surface of the parabolic mirror is curved, so this compensation 
can occur for only one off-axis displacement s. It is suggested that this compensation occurs for the minimum accordion lattice 
period d m \ n = A/2, when the atoms undergo coherent dynamics. The maximum period lattice will still be somewhat distorted; 
however, this is unimportant, as coherent dynamics do not occur in the maximum period lattice. 

Once the accordion lattice has been expanded to its maximum extent the lattice sites should be far enough apart to be resolved 
optically. Rather than directly illuminating the atoms in the accordion lattice, it is suggested that the atoms are at this stage 
loaded into an additional optical lattice — the pinning lattice. This is a deep three-dimensional counter-propagating optical lattice 
(with spatial period A/2); as the intensity of this new lattice is increased, it divides each original accordion lattice well into many 
new wells; atoms are distributed between the new wells (Pig. [3j. The rate of the intensity increase of the pinning lattice should 
be fast enough that there is little transfer of atoms between the pinning lattice sites during the ramp-on of the pinning lattice; 
the distribution of atoms in the pinning lattice sites therefore has roughly the same envelope as the atomic density distribution in 
the accordion lattice before the pinning lattice was raised. After this (non-adiabatic) transfer process, the atoms will be trapped 
at a finite temperature in the pinning lattice. The accordion and pinning lattices do not need to be phase-locked. Once in the 
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FIG. 3: Schematic overview of the accordion and pinning lattice measurement technique. The measurement sequence starts in 

the short-period accordion lattice (1) once the quantum dynamics of interest have taken place within the ultracold atomic 
sample (black). The accordion lattice height is increased, followed by the lattice period (2). The deep pinning lattice is raised 
(3); atoms in the same well in the accordion lattice are predominantly trapped in different wells in the pinning lattice. When the 
sub-Doppler cooling beams are turned on (4), the atoms fluoresce without being lost to light-assisted collisions. For certain 
pinning lattice and cooling beam parameters the mean lifetime of atoms at pinning sites can be many seconds ||4j]. While the 
individual pinning lattice sites cannot be resolved, atoms from different accordion lattice sites should be distinguishable using a 

microscope objective lens. 



pinning lattice, atoms from the same accordion lattice site will be, to a high probability, isolated from each other. When the 
atoms are then probed with resonant light they cannot undergo light-assisted collisions; this addresses the third issue discussed 
in the introduction. Note that the probing lattice can be much closer to the atomic resonance (e.g. a few thousand linewidths [0]) 
than the accordion lattice, which should be far-detuned to minimize light scattering during the coherent dynamics. 

The probability of mutually isolating the atoms depends on the ratio of the lattice site densities (2d max /A) 3 in the expanded 
accordion and pinning lattices; for example this ratio is 10 3 for d max = 5A. Defining n a as the number of atoms at an accordion 
lattice site, and («,) and («,) as the average number of atoms at pinning lattice site i (which overlaps with the original accordion 
lattice site) before and after light-assisted collisions, it is found that 




The fidelity of the division process is defined as the ratio of the average number of atoms observed at an accordion lattice site 
(after light-assisted collisions) to the true (original) number of atoms at that site; it is found to be 

" i 

for («;)«: 1. For example, assuming the atoms are at a temperature of Vo/(4A:b) in an accordion lattice with spatial period 
^max = 5/1 and potential barrier Vo between sites (and that there is negligible transfer of atoms between pinning lattice sites 
during the ramp-on of the pinning lattice), the division fidelity is approximately 98.4% for n a = 2, 96.9% for n a = 3, and 95.3% 
for n a = 4. Increasing the temperature of the atoms in the accordion lattice before loading into the pinning lattice actually 
improves the measurement fidelity, as the atoms are distributed between more pinning lattice sites. 

Once the atoms are in separate sites of the pinning lattice they are subjected to polarization gradient cooling in three di- 
mensions; the scattered light is collected to detect the atoms. The dynamics of this system, together with appropriate optical 
configurations and parameter ranges, are discussed in Ref . J3l ■ There, it is shown that a parameter range exists in which individ- 
ual atoms are confined to separate pinning lattice sites, while scattering photons at rates of around 10 6 s _1 , until background gas 
collisions remove the atoms from the sample. To measure the three-dimensional density distribution of a sample tomographic 
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FIG. 4: Simulation of isolated Rb atoms in the pinning lattice undergoing polarization gradient cooling on the D2 transition 
while being illuminated with probe light on the Di transition. The lifetime is the average time for an atom at a pinning lattice 
site (initially at equilibrium temperature) to leave that site(see Ref. J3l for details). The parameters used for the simulation are 
the same for the ID alternating case in Table 1 of Ref. |4J], with the exception of the frequency and intensity of the pinning 
lattice light, which has a frequency detuning of 2000 linewidths to the blue of the F = 2— > 3' D2 line, and an intensity of 
1.4 x 10 4 mWcirT 2 per beam; this gives the (bare) pinning lattice a depth of 410 [xK. In addition, the simulation includes 
30 mWcirT 2 of repumping light resonant with the F = 1 — > 2' D2 transition, and probing light in a linearly polarized beam 
resonant with the F = 2— » 1' Di (795 nm) transition. The error bars indicate the standard deviations of the quantities; these are 
due to the stochastic nature of the Monte Carlo simulation. The shaded region illustrates the approximate range of Dj 
intensities which would be useful for the measurement technique described in this article. 



measurements are required, with multiple images taken at different depths within the sample 15J. This can be achieved by 
translating the pinning lattice between images using electro-optical or piezoelectric means. 

For one-dimensional, two-dimensional and thin three-dimensional samples, small clusters of atoms corresponding to the 
sites of the original accordion lattice will be resolvable using the above method. Notably, the method offers advantages over the 
method of the experiments described in Refs. J2l and 13] by enabling the use of a smaller period (A/2) lattice, which has a smaller 
site volume and so an enhanced non-linear interaction parameter during the coherent quantum dynamics; and by measurement 
of the true number of atoms at a site, rather than the parity of this number. 

For thicker three-dimensional samples, the fourth issue discussed in the introduction to this article becomes important; if 
the fluorescent scatter from out-of-focus atoms is too great, the signal from atoms in the focal plane of the imaging system 
may be degraded or rendered indistinguishable from the background light in the image. Whether this is the case in a particular 
experiment depends on the density and thickness of the sample, the spatial period of the expanded accordion lattice, and the 
numerical aperture of the imaging system. 

To suppress background light from out-of-focus atoms in thick three-dimensional samples an additional 'probe' beam can be 
used. It is intended that the probing beam illuminates the sample in a thin sheet co-incident with the focal plane of the imaging 
system. The scattered probe light is collected to form the images, rather than the cooling light. Conveniently, the probe beam is 
resonant with a different atomic transition than the cooling and pinning lattice light, so that all light except that scattered from 
the probe beam is efficiently filtered out of the signal. Scattered light from out-of-focus atoms is therefore highly suppressed in 
the resulting images. 

The probe beam changes the polarization gradient cooling dynamics of atoms in the pinning lattice. The effect of the probe 
beam is investigated using the hybrid Monte Carlo-Master Equation technique described in Ref. fl. The simulation technique 
described in the reference is altered for the current situation by the addition of further laser beams and atomic states (with full 
orientational structure), to incorporate probe beam scattering and optical repumping dynamics. 

The simulation results are shown in Fig.|4]as a function of the probe beam intensity. The probing light disrupts the polarization 
gradient cooling process by mixing and changing the relative energies of the ground state orientational sublevels; this is similar 
to the inhibition of polarization gradient cooling in a state-dependent lattice or zero velocity cooling in a magnetic 

field lllh . Consequently the temperature of atoms in the pinning lattice increases, and the lifetime of an atom at a pinning 
lattice site decreases, with increased probe light intensity. If atoms leave their pinning lattice sites during the measurement, the 
measurement is degraded; this sets an upper limit on the useful probe beam intensity. It is seen from Fig. [4] that a probe beam 
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scattering rate of around 10 5 s 1 is possible with a mean site lifetime of order 100 s; this lifetime is more than adequate for a 
probe beam exposure time of around 0. 1 s to Is. 

In order to demonstrate the potential of this technique it is worth illustrating these arguments with a concrete example of 
a possible lattice and probe beam geometry. It is assumed that a three-dimensional parabolic accordion lattice design is used 
with a mirror element similar to that shown in Fig. [2] The focal length of the parabola is assumed to be 8 mm. The minimum 
distance between the accordion lattice beams 2s m i n is assumed to be 1.2 mm (for the four horizontal incoming beams of Fig. [2]) 
and 0.6 mm (for the two vertical beams of Fig. |2j, so that for a laser wavelength of 830 nm the primitive cell of the expanded 
accordion lattice measures (5.5x5.5x11.1) um. It is advantageous to observe the lattice along a direction perpendicular to a 
diagonal lattice plane, as this increases the distance between lattice sites along the line of sight. Here it is assumed that the focal 
plane of the imaging system lies parallel to the (1 12) lattice plane (consequently the line of sight is accessible and lies in the 
horizontal plane). The projected distance in the focal plane between neighbouring accordion lattice sites is 4.5 um and 3.9 um 
in the horizontal and vertical directions respectively. The accordion lattice sites are therefore easily resolved with an imaging 
system which has a resolution of 1.3 um (NA = 0.35 at 795 nm). It is clear that the pinning lattice sites are not resolved; however 
this is not required. The distance between the expanded accordion lattice sites along the line of sight is 19.2 um. The probe 
beam sheet has a waist that is several hundred micrometres in the thick direction and 3.5 [xm in the thin direction; this grows to 
4.9 um one Rayleigh length (48 [xm) along the beam from the beam waist. The probe beam, which is parallel to the (112) lattice 
plane, therefore highly discriminates between different lattice sites along the line of sight. The probe beam alignment should be 
sufficient that the beam remains approximately co-incident with the focal plane of the imaging system across the width of the 
image; in comparison, the focal depth of a NA = 0.35 system is around 5 urn at 795 nm. Many exposures are taken as the pinning 
lattice is translated along the line of sight of the imaging system, with the exposures spaced by around half the probe beam waist 
(around 2 um in this case); in this way the site-by-site three-dimensional distribution of atoms in the original accordion lattice 
can be reconstructed. 

In order to accurately reconstruct the atomic spatial distribution the fluorescence signal from each atom should have known 
amplitude. While predictable variations in the signal per atom, for example due to intensity envelope of the probe beam, may be 
known and compensated, unpredictable variations in the scattering rate should be minimized. 

One source of variability in the probe beam photon scattering rate is caused by differences in the relative phases of the cooling 
beams and the probe beam at the pinning lattice site. As discussed in Ref. 0], changing the relative phases of the cooling beams 
at a lattice site can change both the cooling properties and the photon scattering rate of an atom at that lattice site. This is true 
even when a separate probe beam is used; changing the relative phases of the cooling beams alters the light field at a lattice 
site, which changes the atomic state occupancies and the probe beam scattering rate. For this reason it is suggested that the ID 
alternating cooling beam configuration is used (see description in Ref. J3l) with an additional small frequency offset between 
the counter-propagating beams. The frequency offset rotates the net polarization vector of the (one-dimensional) cooling light 
at the lattice site relative to the probe light polarization; the scattering rate from the probe beam, when averaged over complete 
cycles, is now independent of all the phases of the beams at the lattice site. For the simulation of Fig. |4] a frequency difference 
of 3.6 kHz was used between the counter-propagating beams ( 1/5 of the flash cycle rate), reducing the standard deviation in the 
scattering rate between lattice sites from 9% (without this frequency offset) to negligible levels (less than 1%). 

Another source of variability in the scattered signal is due to the partially coherent nature of the scattered light. The coherent 



fraction of the total scattered light is IU2I1 



2l«^»l 2 /Z <<Z W >); (3) 



here is the atomic raising (lowering) operator, the index e runs over {cr + , <x~, n), and the doubled angled brackets denote time 
averaging of the expectation values. The coherent fraction is found to be 14% for the simulations of Fig.|4] This is lower than for 
a two level atom under monochromatic illumination due to the multiple spontaneous emission routes and due to mixing of the 
atomic population between ground orientational substates by the cooling light. Interference between the coherent components 
of the scattered light from two or more unresolved atoms (i.e. atoms at the same original accordion lattice site) causes anisotropy 
in the scattered radiation and so unpredictable variations in the observed signal. These variations can be avoided by observing 
the atoms along the polarization vector of the probe light; this greatly diminishes the coherent proportion of the scattered light 
(for example, from 14% to 1.3% assuming an objective lens of NA = 0.35). 

An unavoidable source of variation in the observed signal is the photon shot noise. To take an example, 100 photons are 
observed per atom during an exposure time of 0.1 s for a scattering rate of 10 5 s~', assuming 1% of the scattered photons 
contribute to the observed signal; the shot noise is therefore 10% for 1 atom, 6% for 2 atoms, etc. This level of shot noise is 
small enough that to 5 atoms per lattice site can be reliably counted (assuming other noise sources can be suppressed). A 
tomographic measurement of the atomic density of an extended sample may require between 20 to 40 images; the complete 
measurement would take 2 s to 4 s to complete at this frame rate. 
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The reabsorbed fraction of the scattered radiation depends on the optical depth of the sample l/ncr, where n is the atomic 
density and cr is the radiative cross-section. For example, a sample in an accordion lattice (with maximum spatial period 
(5x5x10) (im) with an average of one atom per site, probed with the parameters used in Fig. [4] would have an optical depth of 
8 mm for the probe light and 15 cm for the cooling light. These optical depths are large enough that the effects of reabsorption 
are negligible. 

Although the technique is designed to measure the site occupancies of atoms in an optical lattice, it is also suitable for high 
resolution atom-by-atom measurements of any ultracold atomic sample which can be loaded into the accordion or pinning lattice 
(from, for example, a magnetic trap or free space) prior to the measurement without excessive distortion of the atomic density 
distribution. Site-selective preparation of a two-dimensional atomic ensemble would also be possible in the expanded accordion 
lattice, for example with use of an extra tightly-focussed photoionizing laser beam to eject atoms at particular accordion lattice 
sites. However it is probable that after this preparation process, extra in-lattice cooling would be needed for atoms in the short- 
period accordion lattice, in order to transfer them to the ground vibrational state at each site before coherent quantum tunneling 
dynamics take place. 

In conclusion, it has been discussed how an accordion lattice can have a minimum spatial period of A/2. When expanded, each 
accordion lattice site is optically resolvable. A separate pinning lattice can be used to mutually isolate atoms in the expanded 
accordion lattice, so that atoms cannot be lost due to light-assisted collisions during fluorescence imaging. It has been shown how 
a separate probing beam can be used to aid tomographic measurements of extended three-dimensional samples. Taken together, 
a technique is outlined which is capable of measuring the three-dimensional atom-by-atom site-resolved spatial distribution of 
ultracold atoms in a short spatial period (A/2) optical lattice. 
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